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ABSTRACT 
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V a signal- to -noise  ratio  between  0.2  and  0.1  on  the  original  single-.  ^ ) 
component  records,  amplitudes  for  Rayleigh  waves  over  oceanic  paths  of  155° 

I / '''' "\ 

at  station  MAT  and  98°  at  station  KIP  have  been  determined  as  12  mp  and 
1 24  mU  peak-to-peak,  respectively,  with  a standard  error  of  less  than 
11  per  cent.  In  each  case  the  processed  correlation  signal  is  the  highest 
in  a half-hour  record.  The  method  makes  use  of  preliminary  high-pass 
filtering  and  normalized  reference  earthquake -matched  filtering,  and  takes 
full  advantage  of  the  well-dispersed  oceanic  surface  wave. 

The  method  also  provides  high  resolution  of  co-located  events  with  short 
tinr  separation,  or  of  widely  spaced  events  with  Rayleigh  waves  arriving 
nearly  simultaneously  at  a single  station,  when  the  summed  vertical  and 
radial  matched  filtered  components  are  used. Examples  include:  (1)  clear 

separation  and  amplitude  determination  at  stations  KIP  and  MAT  of  two 
Ms  = 6.5  earthquakes  located  0.7°  and  145  sec  apart  off  the  coast  of  central 
Chile;  (2)  clear  separation  at  station  KIP  of  a Novaya  Zemlya  mb  = 4.8  event 
from  interfering  Rayleigh  waves  of  an  mb  = 5.0  Kermidec  Island  earthquake 
arriving  120  to  140  sec  prior  to  the  searched  event,  with  almost  complete 
elimination  of  interference  on  the  summed  vertical  and  radial  processed 
components;  and  (3)  clear  separation  at  station  KIP  of  two  co-located 
m,  = 4.4  and  4.5  earthquakes  6 min  apart  off  the  coast  of  Chile,  with 

D 

determination  of  their  amplitudes  in  the  presence  of  interfering  Rayleigh 
waves  from  two  central  Alaska  earthquakes,  the  first  (mb  =4.1)  arriving 
15  min  prior  to  the  first  Chile  Rayleigh  wave  and  the  second  between  the  two 

Chile  arrivals. 
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The  single-station  threshold  reached  (10  and  25  digital  units,  p-p)  for 


stations  MAT  and  KIP  at  155°  and  98°,  respectively,  corresponds  to  an 


M = 3.3  and  probably  can  be  improved  further.  Beam  focusing  is  obtained  by 


referring  the  individual  station  correlation  (reference  earthquake  matched 


filter)  outputs  to  the  origin  time  and  summing  these  output  from  the  array 


of  the  randomly  spaced  HGLP  stations. 


It  is  shown  that  the  computer  determined  magnitudes  (relative  to  the 


reference  event)  are  very  stable  among  different  components  and  different 


stations  and  varying  by  less  than  0.1  magnitude. 


Automated  amplitude  and  phase  response  of  the  complete  seismometer- 


recording system  is  obtained  from  step  inputs  to  the  calibration  coil.  High 


accuracy  is  achieved  by  summing  as  many  pulses  as  desired  (to  eliminate 


background  noise)  by  a correlation  technique  and  subsequent  Fourier  analysis. 


The  only  parameters  required  are  the  seismometer  mass,  the  Cal-coil  constant 


(referred  to  the  center  of  mass  if  appropriate)  and  current,  and  the  precise 


onset  time  of  one  reference  calibration  current,  which  are  all  very  stable 


over  long  time  periods.  Application  to  the  High  Gain  Long  Period  system  at 


KIP  yields  the  magnification  curve  from  only  six  pulses  with  less  scatter 


(<  ± 5%  for  periods  larger  than  20  sec)  than  routine  steady-state  calibrations, 


Deconvolution  of  the  digital  seismograms  results  in  retrieval  of  the 


ground  motion  (in  the  frequency  range  of  interest)  by  the  use  of  the  complex 


Fourier  coefficients  obtained  from  the  calibration  method. 


All  computer  programs  developed  or  used  for  purpose  of  this  work  are 


presented  in  a (separate)  report. 
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j I n t r o d uct ion 

A number  of  carefully  designed,  3-coraponent  High  Gain  Long 
Period  (HGLP)  seismometer  stations  have  been  installed  in  sever 
locations.  This  network  is  aimed  at  exploring  the  worldwide,  stable 
earth-noise  background  window  between  30  and  40  sec  period  and  lower- 
ing the  threshold  of  detectable  surface  waves.  Gains  from  ground 
motion  to  record  trace  amplitude  for  "good"  stations  (OGD  and  CTA) 
are  typically  over  100,000  in  the  frequency  band  of  the  noise  window, 
whereas  those  stations  located  in  areas  of  high  oceanic  microseism 
(like  KIP)  have  magnifications  from  20,000  to  40,000.  Many  of  these 
HGLP  stations  are  located  in  and  around  the  Pacific. 

The  reported  noise  background  of  the  better  stations  (15  to 

25  m^  at  40  sec  period)  is  such  that  earthquakes  with  Mg (T=40)  = 2.1 

should  be  detected  at  a distance  of  25°,  those  with  Ms(T=40)  = 3.1 

near  90°.  This  is  based  on  the  40  sec  surface  wave  magnitude  deter- 

1 6 6 

mination  using  Gutenberg's  attenuation  as  A (A  = distance  in 

degrees)  beyond  A = 25°  and  assuming  that  s i gn a 1- to-no is e ratio  on 
the  record  is  near  1:1  for  visual  detection  of  the  surface  waves  by 
a trained  observer.  The  magnitude  of  detectable  earthquakes  ill 
be  up  to  0.4  higher  for  the  HGLP  stations  with  the  lesser  gains. 

In  the  original  proposal  the  Pacific  HGLP  stations  were  consi- 
dered as  a randomly  spaced  array  and  we  had  suggested  a combination 
of  beam  steering  and  signal  extraction  methods  in  the  presence  of 
high,  but  uncorrelated  noise,  in  order  tu  lower  the  detection  thres- 
hold for  small-magnitude  sh al low- f ocu s earthquake  ar  urn  he  Pacific. 

It  was  further  suggested  to  study  the  variation  of  Rayleigh  wave  ra- 
diation pattern  for  different  periods  and  for  the  same  source  areas, 
and  investigate  attenuation  of  40  sec  Rayleigh  waves  across  the  Pacific 
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This  study  was  aimed  at  testing  regional  focal  mechanism 
stability  of  shallow  focus  earthquakes  of  small  magnitudes  so 
that  reference  signals  at  a given  station  in  the  m range  from 
5 to  5.5  could  be  used  for  extracting  signals  of  earthquakes 
of  smaller  magnitudes  by  cross-correlation  methods  that  are 
more  promising  when  the  signals  are  well  dispersed.  Oceanic 
Rayleigh  wave  signals  are  well  dispersed  for  periods  shorter 
than  25  to  28  sec  (01 iver , 1962;  Kuo  et  al. , 1962)  especially 
when  the  generally  long  travel  path  from  c ir cum-Pacif ic  earth- 
quake waves  to  the  KIP  station  are  considered.  The  cross 
correlation,  in  addition,  when  normalized  to  the  reference  event 
signal  will  yield  the  magnitude  for  the  small  quake  and  can  be 
based  on  a selected  frequency  band. 

Other  methods  succeeding  in  s igna 1- t o -no i s e ratio  improvements, 
such  as  polarization  filtering  of  Love  and  Rayleigh  waves  (Simmons , 
1968;  Choy  and  McCamy , 1973)  were  not  considered  at  this  point. 

The  regionalization  approach  avoids  the  uncertainties  invol- 
ving the  distortion  of  the  Rayleigh  waves  spectrum  due  to 
frequency,  path  and  station  site  dependence,  since  only  events 
from  the  same  source  area  at  a given  station  are  compared,  a 
point  recently  stressed  by  Alexander  (1973) 

Preliminary  visual  inspection  of  che  high  gain  photographic 
records  from  KIP  for  the  first  one-half  to  two-thirds  of  each 
of  the  months  of  Spetember,  October  and  November  1973  indicated 
that  many  closely  matched  Rayleigh  wave  trains  could  be  found 
that  would  be  suitable  for  digital  processing.  These  included 
earthquakes  from  off  the  coast  of  Chile,  off  the  coast  of  Mexico, 
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Gulf  of  Mexico,  Gulf  of  Alaska,  Kurile  Island  area,  Luzon- 
phillpplne,  Fiji,  south  of  Kermadec  and  the  mid-Atlantic  idge. 
Preliminary  locations  were  obtained  from  LASA/Nors  ar  bulletins. 
An  example  of  matching  signals  from  the  Fiji  area  is  shown  in 


Fig.  1 for  mB  (LASA)  from  3.9  to  4.5. 


As  a result  a digital  tape  and  computer  plotted  records  for 
selected  events  were  obtained  from  the  Albuquerque  S eismological 
Center  (ASC) , aimed  primarily  at  the  use  for  establishing  and 
verifying  the  computer  programs.  For  further  check  on  matching 
of  Rayleigh  waves,  earthquakes  from  the  south  of  Kermadec  area 
(on  October  15,  1973;  mfi  from  5.1  to  3.7)  had  been  included, 
only  the  larger  shock  was  given  in  the  LASA  bulletin.  The 
smaller  shocks  could  be  identified  by  their  almost  perfect 
match  of  surface  waves  over  10  to  15  minutes.  The  time  differ- 
ence of  occurrence  was  determine^  from  this  match  (±  1 sec)  and 
later  confirmed  from  the  NOAA/USGS  PDE  cards.  Use  of  Everden's 


rr 


M = log  — + 1.00  log  A + 0.92 
s T 

yields  Ms(38)  - 4.4  for  the  largest  and  Mg  - 3.6  for  the  smallest 
shock.  Gutenberg's  formula  results  in  a 0.2  higher  magnitude. 

It  was  therefore  anticipated  that  shocks  of  Mg  - 3.0  to  3.3 
could  be  identified  for  this  area  in  the  normal  background  noise 
of  the  KIP  station.  Similarly  Kurile_J^^  earthquakes  near 
45. 8N  and  151. 7E  could  be  recognized  down  to  Mfi  = 4.3  (4  October 
1973  origin  time  G.S.  185701.6,  Mg  = 4.3,  normal  depth,  p p 
amplitude  on  record  3mm  (Mg  (calculated  from  M£)  - 2.9)  at  a 

distance  near  5200  km).  However,  comparing  the  P-P  amplitude 
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to  the  3 October  1974  quake  (p-p  amplitude  on  photographic 

record  84  mm)  yields  a M of  3.1. 

s 

Figure  3 shows  the  (signal  + noise)  /noise  power  density 
(smoothed  over  5 data  points)  for  the  Kurile  Island  earthquake. 
The  vertical  component  had  instrument  correction  included  and 
is  the  same  as  in  Fig.  2 whereas  the  horizontal  component  did 
not  include  instrumental  corrections.  Signal  to  noise  power 
exceeds  300.  Assuming  signal  to  noise  power  near  1 for  detec- 
tion the  threshold  magnitude  would  be  1.24  lower  than  the  one 
for  this  3 October  1973  Kurile  Island  earthquake  or  near  3.3 
(Mg)  which  is  in  good  agreement  with  visual  results  of  the 
preceding  section. 

No  cr o s s - corr e 1 a t ion  methods  had  been  used  so  far  but  an 

improvement  of  between  0.3  to  0.5  in  M of  the  threshold  was 

s 

expected . 

The  power  spectra  also  suggested  a prefiltering  to  elimi- 
nate noise  at  periods  (longer  than  40  to  50  sec.)  where 
essentially  no  useful  signal  was  found. 
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23  OCT  19/3  FIJI  ISLAND  REGION 


Fiji  Island  earthqaukes  redrawn  from  photographic  records.  Nearly 
identical  signals  indicate  constant  focal  mechanism. None  were 
listed  in  either  the  NORSAR  bulletin  or  the  U.S.G.S.  monthly 
listings.  Four  are  given  by  the  LASA  bulletin. The  smallest  one 
( not  listed  ) has  a surface  wave  magnitude  of  3.?  . 


Power  density  spectrum  for  the  vertical  component  at  KIP  of  the 
Kuril  Island  to  (20)=4.5  earthquake^  Oct.  1973  and  a preceding 
equally  long  noise  section -(l024sec) (left  hand  scale) . Squares 
indicate  (S+ft)/N  ratio  (right  hand  scale) . Ins t rumen!  response 
was  included. 


x EW 


(S+N)/N  power  ratio  for  ail  three  components  of  the  same  earth- 
quaKe  as  in  previous  figure,  horizontals  however  without  instru 
ment  corrections. 


II  Ac  comp li shments 


The  major  results  from  the  research  under  this  grant  have 
been  either  published  (Berg,  1974;  Berg,  1975)  , are  in  press 

(Berg  and  Chesley,  1976;  Chesley  and  Berg,  1976)  or  will  be 
reported  in  the  later  section. 

The  first  two  publications  (Berg,  1974  and  1975)  essentially 
describe  the  extraction  and  amplitude  determination  from  the 
high  gain  long  period  stations  by  the  use  of  matched  reference 
earthquake  filters  applied  at  a single  station.  Therefore, 
only  the  Abstract  of  one  (Berg,  1975)  is  given  (Section  IIA‘. 

The  paper:  Automated  High-precision  Amplitude  and  Phase 

Calibration  of  Seismic  Systems  (Berg  and  Chesley,  1976)  will 
appear  in  the  Bulletin  of  the  Sei smologic al  Society  of  America 
(probably  August  1976) . As  this  paper  is  not  yet  available  it 
is  given  in  its  accepted  version  (Sect.  IIB) . This  method  was 
developed  to  obtain  accurate  amplitude  and  phase  information 
in  view  of  determining  spectral  amplitude  ratios  automatically 
from  the  calibration  pulses  that  are  available  on  each  HGLP 
data  tape  and  thus  making  also  monthly  sinusoidal  calibration 
unnec  e s s ary . 

The  last  reference  (Chesley  and  Berg,  1976)  will  be  publish- 
ed as  a Hawaii  Institute  of  Geophysics  Research  report  and  con- 
tains the  computer  programs  that  have  been  developed  or  used 
under  this  grant  on  the  IBM  370/158.  For  brevity  only  the 
abstract  is  given  (Sect.  IIC). 
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The-  section  on  "Beam  Focusing"  describes  the  use  of  the 
randomly  spaced  HGLT  station  network  (or  others)  as  an  array 


to  e tract  and  determine  amplitudes  of  Rayleigh  waves  from  very 
weak  events,  with  simultaneous  identification  of  the  event  area. 
The  individual  station  matched  filtered  outputs  as  obtained 
by  the  method  of  Berg  (1974,  1975)  are  focused  on  the  desired 
area,  with  an  example  taken  from  Novaya  Zemlya  (Sect.  IID). 
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Bulletin  of  the  Seismolopica!  Society  ol  America.  Vol.  65,  N0.  6,  pp.  1761-1778.  December  1975 


RAVLflGH  WAVES  FROM  HIGH-GAIN  LONG-PERIOD  STATION'S:  SIGNAL 
EXTRACTION,  AMPLITUDE  DETERMINATION.  AND  SEPARATION  OF 
OVERLAPPING  WAVE  TRAINS 

By  Eduard  Biro 


ABSTRACT 

For  a signal-to-noise  ratio  between  0.2  and  0.1  on  the  original  single-component 
records,  amplitudes  for  Rayleigh  waves  over  oceanic  paths  of  155  at  station 
MAT  and  98  at  station  KIP  have  been  determined  as  12  m/j  and  24  m/i  pcak-to- 
peak,  respectively,  with  a standard  error  of  less  than  11  per  cent.  In  each  case  the 
processed  correlation  signal  is  the  highest  in  a half-hour  record.  The  method  makes 
use  of  preliminary  high-pass  filtering  and  normalized  reference  earthquake- 
matched  filtering,  and  takes  full  advantage  of  the  well-dispersed  oceanic  surface 
wave. 

The  method  also  provides  high  resolution  of  co-Iocated  events  with  short  time 
separation,  or  of  widely  spaced  events  with  Rayleigh  waves  arriving  nearly 
simultaneously  at  a single  station,  when  the  summed  vertical  and  radial  matched 
filtered  components  are  used.  Examples  include:  (1)  clear  separation  and  amplitude 
determination  at  stations  KIP  and  MAT  of  two  Ms  = 6.5  earthquakes  located 
0.7  and  145  see  apart  off  the  coast  of  central  Chile;  (2)  clear  separation  at 
station  KIP  of  a Novaya  Zemlya  mh  = 4.8  event  from  interfering  Ravleigh 
waves  of  an  mh  ~ 5.0  Kcrmadec  Island  earthquake  arriving  120  to  140  sec  prior 
to  the  searched  event,  with  almost  complete  elimination  of  interference  on  the 
summed  vertical  and  radial  processed  components;  and  (3)  clear  separation  at 
station  KIP  of  two  co-located  mh  = 4.4  and  4.5  earthquakes  6 min  apart  ofT  the 
coast  of  Chile,  with  determination  of  their  amplitudes  in  the  presence  of  inter- 
fering Rayleigh  waves  from  two  central  Alaska  earthquakes,  the  first  ( mb  ~ 4.1) 
arriving  15  min  prior  to  the  first  Chile  Rayleigh  wave  and  the  second  between 
the  two  Chile  arrivals. 

The  single-station  threshold  reached  (10  and  25  digital  units,  pp)  for  stations 
MAT  and  KIP  at  155  and  98  , respectively,  corresponds  to  an  Ms  = 3.3  and 
probably  can  be  improved  further. 


1 


p 


[. 1 


II  B AUTOMATED  HIGH-PRECISION’  AMPLITUDE 

AND  PHASE  CALIBRATION  OF  SEISMIC  SYSTEMS 
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Eduard  Berg  and  Duncan  M.  Chesley 

i 


ABSTRACT 

Automated  amplitude  and  phase  response  of  the  complete  seismometer- 
recording  system  is  obtained  from  step  inputs  to  the  calibration  coil. 

High  accuracy  is  achieved  by  summing  as  many  pulses  as  desired  (to 
eliminate  background  noise)  by  a correlation  technique  and  subsequent 
Fourier  analysis.  The  only  parameters  required  are  the  seismometer  mass, 
the  Cal-coil  constant  (referred  to  the  center  of  mass  if  appropriate)  and 
current,  and  the  precise  onset  time  of  one  reference  calibration  current, 
which  are  all  very  stable  over  long  time  periods.  Application  to  the  High 
Gain  Long  Period  system  at  KIP  yields  the  magnification  curve  from  only 
six  pulses  with  less  scatter  (<  ± 5%  for  periods  larger  than  20  sec)  than 
routine  steady-state  calibrations. 
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INTRODUCTION 

Many  papers  have  described  the  calibration  of  seismometer  systems. 

Three  approaches  have  been  developed:  The  parameter  method,  in  which  the 

parameters  necessary  for  calculation  of  amplitude  and  phase  response  have 
been  determined;  the  steady-state  method,  in  which  the  amplitude  and  phase 
response  is  determined  for  discrete  frequencies  either  by  exiting  the 
seismometer  on  a shake  table  or  by  the  use  of  a calibration  coil,  and 
finally  the  transient  method,  in  which  amplitude  and  phase  response  from 
the  recorded  pulse  on  the  seismogram  are  extracted  using  a step  or  impulse 
input  into  a calibration  coil.  Accurate  knowlege  of  the  response  is 
especially  necessary  for  research  on  surface  wave  attenuation  and  dispersion. 
Most  of  the  more  recent  approaches  necessitate  tedious  manual  comparison 
with  pre-calculated  responses  of  calibration  pulses  (Espinosa  et  al. , 1962, 
1965)  or  very  accurate  measurements  of  phase  and  amplitudes  (inducing 
their  temperature  variation)  in  a somewhat  modified  steady-state  response 
approach  (Mitronovas  and  Wielandt,  1975).  Mitchell  and  Landisman  (1969) 
also  developed  a least -square  method  to  determine  free  period  and  damping 
constants  of  an  electromagnetic  seismograph  from  its  calibration  pulse,  by 
weighting  those  frequencies  that  are  the  freest  of  noise  and  calculating 
the  response  for  higher  frequencies  from  the  theory  of  seismometer- 
galvanometer  systems.  It  seems  difficult  however  to  follow  such  an 
approach  when  more  ^aplicated  systems  are  involved.  Such  systems  might 
include  band-rejection  filters,  additional  electronic  filters,  and  digital 
or  FM  long-distance  telemetry  to  the  final  recording  site. 

In  this  paper  the  transient  method  is  used  and  accurate  amplitude 
and  phase  calibration  can  be  most  easily  obtained  where  calibration  pulses 


recorded  on  a routine  (day  to  day)  basis  and  recording  is  on  magnetic 
tape.  By  determination  of  the  exact  start  of  each  individual  recorded 
pulse  by  a correlation  technique  (Berg,  1974  and  1975),  as  many  pulses  as 
desired  for  accuracy  can  be  summed  and  this  sum  Fourier  analyzed.  The 
parameters  necessary  to  obtain  the  complete  absolute  response  are  the 
precise  start  time  of  one  pulse  (a  "ref erence") , the  mass,  the  calibration 
coil  constant  (referred  to  the  center  of  mass  if  appropriate)  and  the 
applied  calibration  coil  current  of  the  seismometer. 

It  should  be  stressed  that  this  method  is  directly  applicable  to  the 
calibration  of  all  systems  equipped  with  a calibration  coil  and  would  be 
easier  than  the  more  cumbersome  steady-state  method.  The  method  also  can 
be  used  (with  the  appropriate  slight  modifications)  in  case  of  capacitor 
plate  calibration  devices  (as  used  in  tidal  gravimeters)  or  weight  lifts 
(as  in  the  Askania-type  bore  hole  tiltmeters)  or  expansion  cells  (and 
pieco  electric  crystals)  as  applied  to  the  Verbandert-Melchior  quartz  tidal 
pendulums.  The  method  is  especially  useful  for  remotely  operating  and 
recording  systems,  such  as  those  in  ocean-bottom,  microearthquake,  or 
planetary  systems,  since  the  parameters  required  (e.g.  mass,  calibration 
coil  constant,  and  current)  are  very  stable. 

It  also  should  be  pointed  out  that  the  total  system  response  is 
obtained,  that  is,  from  the  ground  motion  to  the  final  record.  Such 
systems  might  include  different  transducers,  electronic  amplifiers, 
telemetry  sections  and  filters  before  the  final  record. 

Since  transient  calibration  in  operational  installations  is  performed 
remotely,  no  physical  disturbance  of  the  seismographs  is  required 


In  practice  a ten-day  tape  from  a Long  Period  High  Gain  station 
yields  a sufficient  number  of  pulses  to  obtain  the  response  within 
better  than  5%  over  the  useful  period  range.  Since  the  complex  Fourier 
coefficients  can  be  stored  in  the  computer,  a recorded  event  (to  be 
analyzed)  can  be  directly  convolved  to  obtain  ground  motion. 


THEORY 


A seismogram  can  be  scanned  and  the  start  time  and  amplitude  (with 
standard  error)  of  the  scanned  pulse  obtained  by  a correlation  method 
(Berg,  1974,  1975),  if  a reference  pulse  is  available  with  known  start 
time.  It  is  assumed  that  (as  is  customary  in  routine  operation  of 
seismograph  stations)  all  calibration  pulses  are  generated  by  the  same 
calibration  current  into  the  coil.  Beginning  with  the  start  times,  all 
individual  pulses  are  summed,  and  the  resulting  amplitude  divided  by  the 
number  of  pulses  used,  to  obtain  a representative  calibration  pulse,  thus 
eliminating  noise  by  an  amount  that  depends  on  the  number  of  pulses  used. 
This  improved  the  high  frequency  calibration  in  long-period  systems 
espec ially. 

The  equation  of  motion  of  the  seismometer  mass  is  given  by  (see 
Jarosch  and  Curtis,  1973  or  Gamburzew,  1964) 


<^I  + 2 + *--£? 


where 


X = displacement  of  mass  with  respect  to  its  frame 
T0  = 2 rr/wo  free  period 


r 


eQ  = h coq  damping  coefficient 


Y = ground  displacement  (positive  in  the  same  direction  as  X) 


Consider  a calibration  current  applied  at  t = 0 causing  a step  function 


in  acceleration  (equivalent) 


d^Y  GI 


t " m • H(t) 


where 


G = calibration  coil  motor  constant  (referred  to  the  center  of 
mass  if  applicable) 

I = applied  calibration  current 
M = seismometer  mass 
H(t)  = Heaviside  step  function. 

In  most  short  period  seismometers,  the  calibration  coil  is  wound  uni- 
ax  ially  with  the  transducer  coll  and  equation  (2)  applies  directly.  In  the 
case,  where  the  calibration  coil  is  physically  located  on  the  boom  but  away 
from  the  center  of  mass,  (such  as  for  most  long  period  instruments)  G has 
to  be  referred  to  the  center  of  mass  position  and  is  given  by 


G = G„ 


where 


= calibration  coil  constant 


= hinge  to  coil  distance 
= hinge  to  center  of  mass  distance, 


r 
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Assuming  no  initial  velocity  or  displacement  of  the  mass  and  operating 
on  the  differential  equation  with  the  Laplace  transform  (after  introduc- 
tion of  the  acceleration  term),  '■ 

7 9 GT 

(S  + 2 Cq  S + u)o  } XA  = - - . (3) 

This  equation  reduces  to 

F (icu)2  + 2e  i(u  + 2 ] x (oj)  = - (4) 

L o o J Mxtu 

where  x (w)  is  the  Fourier  transform  of  the  mass  response  to  the  step 
input  in  acceleration. 

For  a sinusoidal  ground  displacement  of  amplitude  A,  the  differential 
equation  yields  the  mass  response  Xq(w) 

[ (i<u)2  + 2eoi0)  + U)o2  ] XD(W)  = - (im)2  A (5) 

and  dividing  the  last  equation  by  the  acceleration  response  equation 
yields 

XD(0J)  = XA(UJ)  ‘ (iOJ)3  ~ (6) 

If  XR  is  the  frequency  response  of  the  recording  equipment  including 
that  of  the  transducer,  the  output  response  to  the  step  input  in 
acceleration  will  be  Fa(<d)  = XA(w)  • XR(w)  and  the  response  to  the 
ground  displacement  will  be  FD(cu)  = XD(u>)  • xR(w)  • 

Since  Fa(uj)  can  be  obtained  through  Fourier  analysis  of  the  recorded 


output  signal  from  the  step  input  in  acceleration,  the  total  system 
response  Fjj(u;)  to  ground  displacement  of  amplitude  A 


can  be  calculated,  provided  M,  G (or  GQ,  r£  and  rcm) , and  I are 
measured.  This  is  true  with  or  without  back  coupling.  This  can  also 
be  seen  by  taking  the  approach  of  communication  theory,  and  considering 
the  seismometer-recording  system  as  a "black  box." 

APPLICATION 

The  method  was  applied  to  obtain  the  calibration  in  amplitude  and 
phase  for  the  Kipapa,  Hawaii  (KIP)  High  Gain  Long  Period  system.  A 
prototype  of  such  a system  was  described  by  Pomeroy  et  ’1.  (1969)  and 
operational  systems  and  installation  techniques  are  dt  ■ ibed  in  a 
Lament -Doherty  Geological  Observatory  (1971)  technical  report.  The  KIP 
system  (in  1974)  includes  a long  period  seismometer  (TQ  « 30  sec)  sheltered 
in  a pressure  tank,  with  a coil  transducer  for  the  seismic  output,  coupled 
to  a galvanometer  (Tg  ~ 100  sec)  in  a photo  tube  amplifier,  and  additional 
filters.  Ttv  aim  of  part,  of  the  filters  is  to  reject  the  predominant 
ocean-generated  6-sec  microseisms  that  strongly  limit  the  magnification 
of  the  WWSN  long  period  components.  The  output  is  then  recorded 
digitally.  The  reference  pulse  was  obtained  by  summing  the  six  pulses  from 
Nov.  22,  24,  and  26,  1974.  Negative -going  pulses  were  added  with  the 
inverse  sign.  The  time  coincidence  was  obtained  from  the  start  time  as 
indicated  by  the  station  log  sheet  and  verified  by  the  correlation  method 
(see  later)  . Figure  1 shows  the  average  pulses  thus  obtained  for  the 
Z (CH:Sl),  NS  (CH:S2),  and  EW  (CH:S3)  components,  where  the  start  time 
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is  at  the  minute  mark  labeled  "2"  and  the  y exis  is  labeled  in  digital 
units.  Start  times  of  the  individual  pulses  were  verified  by  the  correla- 
tion method  of  Berg  (1974,  1975)  and  are  shown  for  the  22  Nov.  pulses  in 
Figure  2.  In  this  figure  the  correlation  traces  are  labeled  CH:C__  and  the 
time  on  the  correlation  traces  is  in  terms  of  start  time  of  the  reference 
traces,  so  that  maximum  correlation  (±  1)  is  at  the  start  time  of  the 
calibration  pulse  in  the  record  section.  Table  1 shows  the  computer 
output  data  for  the  calibration  pulse  correlations  for  Nov.  22,  24,  and 
26,  1974.  The  time  column  gives  the  calculated  start  times  of  the 
calibration  current;  "data  for  maximum"  refers  to  the  positive-going 
pulse,  and  "data  for  minimum"  to  the  negative -going  pulse.  All  start 
times  thus  obtained  from  the  correlation  method  are  within  ± 1 second 
(the  digitizing  interval)  from  those  in  the  station  log  sheet.  The 
column  "slope"  gives  the  relative  amplitude  with  respect  to  the  reference 
(the  average  of  the  six  pulses  for  the  particular  component)  "error"  is 
one  standard  deviation  of  the  "slope"  (relative  amplitude)  and  "amplitude" 
is  the  peak  to  trough  amplitude  in  digital  units.  It  can  be  seen  from 
the  table  that  with  the  exception  of  day  326,  all  positive-going  pulses 
have  larger  amplitudes  than  the  negative -going  pulses,  possibly  pointing 
towards  some  instrumental  asymmetry.  This  is  also  reflected  in  the  fact 
that  the  amplitude  errors  (one  standard  deviation)  are  considerably 
smaller  than  the  actual  deviations  from  the  average  amplitudes.  From 
statistics,  one  would  expect  a 99%  confidence  interval  for  about  2.6 
standard  deviations  (see  Crow  et  al.,  1960,  p.  160  and  231).  The  actual 
average  peak  to  trough  amplitudes  (three  positive-  and  three  negative- 
going pulses)  are:  for  Z,  1100.0;  for  N-S,  1694.7;  and  for  E-W,  1922.0 
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digital  units,  whereas  the  average  amplitudes  of  the  three  pos it  i/e -going 
pulses  (as  obtained  from  the  correlation  output  in  Table  1)  are:  for  Z, 

1120.1  db  8.8;  for  N-S,  1731.2  ± 69.4;  and  for  E-W,  1948.4  ± 77.3. 

The  average  pulses  from  the  three  components  have  been  submitted  to 
a fast  Fourier  transform  to  obtain  the  response  to  a step  of  acceleration 
Fa(W)  (in  equation  7)  and  (by  multiplying  with  (ict)  A*M/(G*I))  the 
displacement  response  Fd(oj).  If  the  record  amplitudes  are  measured  in 
digital  units  (Dig),  CD  in  1/sec,  A in  meters,  M in  kg,  G in  Newtons/Amp 
(mkg/(sec2Amp))  and  I in  Amp,  the  output  is  obtained  in  Dig/m. 

The  instrument  constants  applied  here  were  supplied  by  John  Hoffman 


(Albuquerque  Seismological  Laboratory,  personal  communication)  and  are  as 
follows : 


Mass 

(kg) 

Go 

(Newton/Amp) 

rc 

(mm) 

rcM 

(mm) 

G (calculated) 
“ Go*rc^rcM 

Z 

13.2 

0.0309 

378 

230 

0.05078 

N-S 

10 

0.0237 

378 

227 

0.03947 

E-W 

10 

0.0277 

378 

227 

0.04613 

<rc 

= coil  distance  to  hinge,  rcl^  - center  of 

mass  distant  to  hinge) 

and 

the 

calioration 

current  in  a 

11  cases 

is  constant  at  5 micro  Amp. 

The 

amplitude 

response  thus 

obtained 

for  the 

three  components  of 

the 

KIP 

High  Gain  Long  Period  instruments 

is  shown  in  Figure  3.  In  this 

figure  the  circles  represent  amplitude  response  for  a 256-sec-long  section 


starting  at  the  onset  of  the  calibration  current,  the  small  crosses 
represent  amplitude  response  for  a 512-sec-long  section  starting  at  the 
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same  time  as  the  256-sec  section.  It  can  be  seen  that  the  longer  section 
introduces  considerably  more  ncise  at  shorter  periods.  Also  shown  for 
comparison  are  the  steady  state  calibration  values  at  100,  70,  40,  and 
30  sec,  obtained  during  the  two  immediately  preceding  routine  checks  on 
10  Oct.  and  8 Nov.,  1974  (John  Hoffman,  personal  communication).  It  can 
be  seen  that  the  .alysis  of  the  transient  response  of  six  averaged 
pulses  results  in  an  amplitude  response  within  5%  of  the  smooth  line  in 
the  worst  case  for  periods  longer  than  20  seconds  on  the  Z and  EW 
components  and  longer  than  24  sec  for  the  NS  component  for  a 256 -sec 
record  length,  whereas  the  steady-state  calibrations  at  30  sec  show 
considerably  more  variation  (typically  10%  or  more) . It  should  be  pointed 
out  that  the  variation  obtained  in  the  transient  analysis  can  easily  be 
reduced  by  application  of  more  calibration  pulses,  especially  during  days 
wnen  microseismic  background  noise  is  low. 

The  phase  response  of  the  three  components  is  given  in  Figure  4. 

The  response  is  smooth  to  the  same  short-period  range  as  the  amplitude 
response.  Since  the  time  of  the  calibration  pulses  used  here  the  system 
has  been  changed  from  photo— tube  amplifiers  to  solid  state,  with  the 
solid  state  response  to  be  as  closely  identical  as  possible  to  the  older 
system.  Jon  Peterson  (written  communication,  December  8,  1975)  supplied 
a phase  response  that  matches  (for  the  N-S  component)  ours  'perfectly 
except  for  a factor  of  tt  (or  180°) . The  inverse  Fourier  transform  of 
our  amplitude  and  phase  response  which  is  the  time  response  to  an  impulse 
in  ground  motion  in  the  positive  direction,  yields  the  first  swing  in  the 
positive  direction.  Therefore  the  phase  response  presented  here  should 
be  the  correct  one. 
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DISCUSSION 

The  limitations  of  the  transient  calibration  are  given  by  the  ratio 
of  the  cal -pulse  amplitude  to  seismic  background  noise  and  by  the  accuracy 
in  the  determination  of  M,  G (or  G0,  rc,  rcM) , and  I.  The  percentage 
errors  in  M,  G,  and  I appear  linearly  in  the  absolute  amplitude  results 
and  can  be  made  small,  and  the  change  with  time  in  a typical  installation 
should  be  very  small  (<  0,17o/y).  They  do  not  affect  phase  or  relative 
amplitude.  The  major  error  source  then  is  provided  by  the  seismic 
background  noise,  since  the  useful  signal  amplitude  is  fixed.  As 
demonstrated  here,  the  noise  background  can  be  reduced  by  averaging  a 
number  of  pulses  by  a correlation  technique.  The  advantage  of  the 
method  is  clear:  the  calibration  can  be  performed  on  a routine  basis 

automatically  without  physically  disturbing  the  seismometer  or  the  record- 
ing equipment  and  no  other  human  error  sources  other  than  in  the  measure- 
ment of  M,  G,  and  I will  occur.  It  is  especially  advantageous  in  remote 
installations  where  telemetry  links  are  involved  since  pulses  arc  much 
more  easily  generated  ard  amplitude-controlled  than  are  the  sinusoidal 
signals  used  in  steady  state  calibration,  and  can  be  applied  automatically 
at  regular  time  intervals . The  recorded  tapes  of  remote  stations  can  be 
processed  in  r.  central  location,  thus  avoiding  displacement  of  personnel 
and  equipment.  Since  for  research  the  calibration  can  be  obtained  usually 
from  the  same  data  tape  as  the  seismic  events,  it  is  possible  to  convert 
the  trace  amplitudes  and  phases  directly  to  ground  motion  by  Fourier 
transforming  the  record,  dividing  by  the  response  and  inverse  Fourier 
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transform  back  to  ground  motion.  Needless  to  say  the  method  can  also  be 
applied  to  short-period  and  other  seismic  systems,  recording  gravimeter, 
and  pendulum-type  tiltmeters. 
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Figure  1.  Average  of  the  three  positive-going  and  the  three  inverted 

negative-going  calibration  pulses  of  station  KIP  for  Nov.  22, 

24,  and  26,  197A;  upper  trace,  vertical;  middle  trace,  N-S; 
bottom  trace,  E-U.  Amplitudes  are  in  digital  units  (y-axis) 
and  time  is  in  minutes  (x-axis) . 

Figure  2.  Lines  3 . 6 , and  9 are  the  average  of  six  calibration  pulses  of 
the  vertical,  N-S,  and  E-W  components  (the  reference  pulses). 
Lines  2,  5,  and  8 are  the  seismograms  for  Nov.  22,  1974. 
Amplitudes  are  in  digital  units,  time  ticks  are  2 minutes  apart. 

* 

Lines  1,  4,  and  7 are  the  correlation  between  the  reference 
pulse  and  the  seismogram  for  the  vertical,  N-S,  and  E-W 
components.  Amplitudes  are  the  correlation  coefficients  (from 
+1  to  -1) . The  maximum  and  minimum  of  these  traces  occur  at 
the  time  of  calibration-pulse  onset  (given  in  Table  1). 

Figure  3.  Amplitude  response  of  the  KIP  High  Gain  Long  Period  system 

obtained  by  Fourier  transform  of  the  average  of  six  calibration 
pulses  for  f;  vertical  (3A) , the  N-S  (3B) , and  E-W  (3C) 
components.  Note  the  scatter  of  the  two  steady  state  calibra- 
tions at  30-sec  period  immediately  preceding  the  time  of  the 
present  transient  calibration.  Reliable  calibration  (with 
less  than  51  amplitude  scatter)  is  obtained  for  a 256-sec-long 
record  section  for  periods  longer  than  20  sec  for  the  Z and  E-W 
component,  and  for  periods  longer  than  24  sec  for  the  N-S 
component . 
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Figure  4.  Phase  re;ponse  of  the  KIP  High  Gain  Long  Period  System  obtained 
by  Fourier  transform  of  the  average  of  six  calibration  pulses 
for  the  (4A)  verticalj  the  (4B)  N-S,  and  the  (4C)  E-W 
components. 
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Computer  Programs  for  Seismology: 
Special  Applications  to  the 
High-Gain  Long  Period  Seismic  Network 

by 

Duncan  M.  Chesley  and  Eduard  Berg 


/ bstract : 

This  report  describes  several  IBM  370/158  compu  er  programs  that  have 
been  applied  to  long  period  seismic  data  on  digital  tapes.  Included  are 
reading  routines  for  tapes  generated  by  the  High-Gain  Long  Period  (HGLP) 
stations,  rotation  of  the  horizontal  components,  low-,  band-  and  high-pass 
digital  filters,  correlation  (or  matched  filtering),  summation  of  correlations, 
beam  focusing,  and  ground  motion  retrieval. 

Results  of  applications,  described  else-  teve  or  in  preparation,  include 
extraction  and  amplitude  determination  of  Rayleigh  waves  with  signal-to-noisc 
ratios  near  1 to  10,  precision  amplitude  and  phase  calibration  of  entire 
seismic  systems  from  transducer  to  final  record,  beam  forming  of  matched 
filtered  outputs  froir  the  randomly  spaced  HGLP  station  array,  and  separation 
of  co-located  multiple  events  with  time  spacing  as  short  as  150  sec  or 
separation  of  events  widely  spaced  geographically  but  arriving  nearly 
simultaneously  at  a given  station. 


Report:  HIG  76-3  estimated  75  pages,  March  1976, Hawaii  Institute  of  Geophysics, 

Univ.  of  Hawaii,  in  press. 


Beam  Focusing  of  the  High-Gain  Long-Period  Stations 
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The  method  (Sect  II  A)  of  extracting  well  dispersed  Rayleigh  waves  with 
amplitudes  far  down  in  the  noise  level  worked  well  at  individual  stations 
and  especially  when  applied  to  wave  trains  predominantly  travelling  a long 
oceanic  path.  It  also  was  argued  (Berg,  1975)  that  slight  differences  in 
earthquake  focal  mechanism,  focal  depth  or  location  between  the  reference 
earthquake  and  the  extracted  one  would  deteriorate  the  correlation,  since  the 
frequency-time  functions  of  the  received  signals  are  somewhat  different 
(besides  the  amplitude  scaling  factor).  A corollary  increase  in  the  standard 
deviation  of  the  amplitude  determination  will  result. 

It  was  thought  then,  that  using  a number  of  n stations  and  adding  their 
correlation  signals  would  result  in  an  increase  of  the  signal  to  noise  ratio 
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by  a factor  of  \ n . It  was  found  that  this  is  nearly  correct,  however 
usually  there  seem  to  be  one  or  two  stations  that  show  exceptionally  clear 
correlation  whereas  the  others  are  nearly  at  the  same  time,  but  their 
amplitudes  are  those  of  the  noise  background.  As  an  example  we  have  chosen 
the  27  Nov  1973  Novaya  Zomiya  explosion  as  a reference,  to  scan  the  records 
of  several  other  stations  for  aftershocks  or  cavity  collapse.  Table  1 gives 
the  source  data. 

Figure  1 shows  the  azimuthal  distribution  and  distances  (in  km)  of  some 
HGLP  stations  with  respect  to  the  reference  event  of  27  November  1973.  In 
addition,  the  great  circle  path  directions  of  arrival  for  the  waves  are 
indicated  by  the  short  arrows  near  each  station  abbreviation.  In  Figure  2, 
the  lower  record  sections  are  the  average  matched  filtered  output  from 
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individual  stations  (marked  CH:SM)  given  in  Figure  3 and  the  uppermost  trace, 
is  the  average  of  the  lower  ones  for  the  time  interval  common  to  all  lower 
traces.  Note  the  different  start  times  of  individual  station  traces.  Since 
the  timing  of  the  matched  filter  output  has  been  referred  to  the  origin  time 
(rather  than  the  start  of  the  Rayleigh  wave)  of  the  reference  event,  the 
correlation  for  the  extracted  event  appears  at  its  own  origin  time  (or  slightly 
offset  if  the  epicenter  is  not  the  same  as  for  the  reference).  The  origin 
time  for  the  reference  event  that  was  used  here  was  07:00:00,  the  average 
between  the  NOAA  time  and  that  in  the  NORSAR  bulletin.  The  beam  minimum  is 
at  09:13:50  in  good  agreement  with  the  NOAA  value  of  09:13:51*  (see  table  1). 

The  beam  output  has  about  twice  the  amplitude  of  the  "background."  We 
therefore  expect  that  it  is  possible  to  further  lower  the  detection  (and 
proper  location)  threshold,  and  still  be  able  to  determine  the  average 
surface  wave  amplitude  as  well.  Note  also  that  the  correlation  "noise"  for 
individual  stations  is  near  0.2  whereas  the  beam  shows  the  reduction  by  the 
square  root  of  the  number  of  stations  used.  The  left  side  of  Figure  2 shows 
the  beam  for  only  those  stations,  where  the  individual  vertical,  the  radial 
or  both  components  had  detected  the  presence  and  determined  the  amplitude 
of  the  surface  (Rayleigh)  waves.  The  right  side  has  the  station  CHG  added 
where  detection  was  not  achieved  (or  was  somewhat  early).  Obviously,  if  an 
event  was  detected  by  short  period  seismometers,  one  is  in  the  favorable 
position  to  choose  only  those  LP  stations  that  have  "detected"  the  surface 
waves  for  amplitude  determination,  and  then  for  the  focused  beam  to  confirm 
the  location.  In  our  example  neither  ALQ  (Fig.  3E)  nor  CHG  (Fig.  3D)  had 
given  a single  component  correlation  that  would  have  been  useful  for 
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amplitude  determination.  Figures  3A  through  3C  give  the  individual  station 
outputs  for  KIP,  TLO  and  MAT,  where  at  least  one  component  succeeded  in 
extracting  the  signal.  Note  the  clear  detection  at  KIP,  Figure  3A  and  on 
the  TLO  vert  ical  (second  trace,  Fig.  3B)  and  the  MAT  N-S  (fifth  trace,  Fig.  3C) . 
Table  2 summarizes  the  results  of  the  detection  at  individual  stations,  for 
the  beam  focusing  for  the  stations  KIP,  TLO,  MAT  and  the  beam  focusing  for 
the  stations' KIP,  TLO,  MAT,  and  CHG.  The  origin  times,  correlation  coeffi- 
cients, p-p  amplitudes  and  magnitude  differences  are  those  obtained  by  the 
computer  printout. 

As  the  table  shows  the  average  surface  wave  magnitude  difference  for 
the  Novaya  Zemlya  event  is  2.29  below  the  reference  of  Ms  = 5.5  resulting 
in  an  Mc  = 3.2  and  all  determinations  are  within  ±0.1  magnitude  of  that 
value,  despite  the  presence  of  interfering  Rayleigh  waves  from  an  mb  = 5.0 
earthquake  S of  Kermadec  (see  Berg,  1975,  Fig.  10).  Note  that  the  signal 
levels  are  lower  than  the  noise  level  for  this  extracted  and  amplitude 
determined  event  (see  Figs.  3A-D)  and  correspond  roughly  to  the  50  to  604 
(all  event)  detection  level  of  the  three  stations  as  determined  by  Lambert, 
Brahl  and  Straus  (November  1973).  The  clearly  negative  correlation  of  KIP 
and  the  beam  outputs  in  Figure  3A  and  Figure  2 together  with  the  position 
(M  = 3.2,  mb  = 4.6  LASA  and  BMO)  on  the  Mg  - mb  plot  for  the  area  presented 
by  Savino  et_  aJL  (1971,  p.  8014)  identify  this  event  as  earthquake  or 
collapse . 

All  seismic  traces  in  Figure  3 had  been  high  pass  filtered  to  eliminate 
the  long  period  portion  of  the  spectrum  where  noise  is  usually  high  and  no 
significant  signal  is  expected  for  very  small  events  (see  also  Section  I). 
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The  question  arises  why  detection  has  not  been  achieved  at  some  of  the 
stations,  notably  CHG,  which  is  at  the  shortest  distance  from  the  epicenter. 

Part  of  the  answer  is  probably  found  in  the  complicated  spectra  of  the 
reference  event,  as  received  at  the  stations.  These  amplitude,  spectra  have 
been  determined  for  the  first  512  sec  of  the  reference  event  without  prefiltering 
or  corrections  for  instrument  response  and  normalized  to  "one"  at  42.7  sec. 

These  spectral  amplitude  ratios  are  shown  in  Figures  4A  through  4E  for  the 

stations  KIP,  TLO,  MAT,  CHG  and  ALQ.  The  high  frequency  cutoff  is  due  in 

part  to  the  cutoff  in  time  of  the  lower  velocity  (higher  frequency)  portion 
of  the  dispersed  Rayleigh  wave  train.  It  seems  as  if  detection  and  amplitude 
determination  is  most  easily  achieved  where  the  spectral  amplitude  ratio 
A(T)/A(42)  is  relatively  high  (above  1.6)  and  relatively  smooth  (as  a 
function  of  period)  in  the  most  useful  period  range.  These  cases  also  seem 
to  correspond  to  paths  that  are  the  least  complicated  in  terms  cf  lithospheric 
structure.  Perhaps  detection  and  amplitude  determination  could  be  achieved 
for  the  event  at  CHG  and  ALQ  by  first  band  pass  filter  the  records  around 
70  and  35  sec  £ CHG  (Fig.  4D)  and  ALQ.  The  rather  strong  long  period  parts 
of  the  station  CHG  and  to  some  extend  ALQ  thus  seem  to  confirm  the  presence 

of  a continental  type  wave  guide  that  has  been  used  by  McDonald  e£  a_l.  (1974) 

for  detection  of  weak  central  Asian  events. 
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Upper  trace  beam  output,  the  average  sum  of  the  lower 
traces,  which  are  the  individual  station  SM  traces  of 
Fig  3.  The  negative  correlation  (with  respect  to  the 
explosion  waves  used  as  reference)  indicates  that  the 
event  was  not  an  explosion  and  gives  the  original  time 
as  09:13:50  and  the  location  near  the  reference  event. 
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KIP  - trace  1:  average  sum  (SM)  of  correlation  traces, 

note  negative  peak  at  09:13:51;  trace  2:  5 matched 

filter  outputs  (Cl,  C2)  for  the  vertical  and  N-S  (radial) 
components;  trace  3:  6 scanned  data  after  HP  prefiltering 

at  45  sec  (Al,  A2)  for  the  vertical,  N-S  component; 
trace  4:  7 reference  explosion  Rayleigh  waves  prefiltered 

(UP  45  sec)  (Al,  A2)  for  the  vertical  and  N-S  (radial) 
components.  For  values  of  correlation,  exact  time  and 
amplitude  determination  see  Table  2. 
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Fig.  4B.  TLO : Vertical  and  Radial  Spectral  Amplitude  ratios 

normalized  at  42.7  sec.,  Novaya  Zemlya  Start  Time 

7:13:00,  Data  length  512  sec. 
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Fig.  4c . MAT:  Vertical  and  N-S  Spectral  Amplitude  ratios 

normalized  at  42.7  sec.,  Novaya  Zemlya 
Start  time  7:19:45,  Data  length  517  sec. 
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Fig.  4D . GHG : Vertical  and  N-S  Spectral  Amplitude 

ratios,  normalized  at  42.7  sec,  Nov ay a Zemlya 
Start  time  7:19:45,  Data  length  512  sec. 
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ratios,  normalized  at  42.7  sec.,  Novaya  Zemlya 
Start  time  7:29:20,  Data  length  512  sec. 
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